Colloidal charge reversal: dependence on the ionic size and the electrolyte concentration 
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Extensive Monte Carlo simulations and scaling arguments are used to study the colloidal charge reversal. 
The critical colloidal surface charge density CTc at which the reversal first appears is found to depend strongly 
on the ionic size. We find that Oc has an inflection point as a function of the electrolyte concentration. The 
width of the plateau region in the vicinity of the inflection point depends on the temperature and the ionic radius 
a. In agreement with the theoretical predictions it is found that the critical colloidal charge above which the 
electrophoretic mobility becomes reversed diverges as Zc ~ in the limit a ^ 0. 



I. INTRODUCTION 

A common way to stabilize colloidal suspensions against 
flocculation and precipitation is by synthesizing particles 
with acidic groups on their surface. When placed in wa- 
ter, these groups dissociate and colloids acquire a net nega- 
tive charge. In aqueous suspension containing only monova- 
lent electrolyte, the long range Coulomb repulsion prevents 
the colloids from approaching one another to distances for 
which the short range van der Waals interaction can lead to 
an iiTeversible sticking and precipitation. However, when be- 
sides the 1:1 electrolyte, suspensions contain some multiva- 
lent counterions, a number of curious and very counterintu- 
itive effects can take place | IJ. For example, it has been ob- 
served that in such suspensions two like-ch^ged colloidal par- 
ticles can attract one another H H B H H, HE H [H lil 
This attraction is not a result of the van der 
Waals interaction, but rather a consequence of strong posi- 
tional correlations between the multivalent counterions sur- 
rounding the colloidal particles |17, 18, 19, 20, 21]. The 
mechanism of this attraction has been studied extensively, yet 
a fully predictive theory of this phenomenon still remains elu- 
sive. 

Another curious effect observed in dilute colloidal suspen- 
sions containing multivalent counterions is the reversal of the 
electrophoretic mobility 1 1, 22, 23, 24]. Since the bare charge 
of colloidal particles is negative, when the electrostatic po- 
tential gradient is established in the suspension, one natu- 
rally expects that the particles should move in the direction 
opposite to the established electric field. Yet, what is of- 
ten found is quiet the opposite — particles drift in the di- 
rection of the field lf25l l26l l27ll . The reversal of the elec- 
trophoretic mobility is a consequence of strong electrostatic 
interaction between the colloidal particles and the multiva- 
lent counterions 1281 i29i l30ll . As a consequence of this cou- 
pling, some counterions become associated (condensed) with 
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the colloidal particle. The positional correlations induced by 
the electrostatic repulsion between the condensed counteri- 
ons can lead to colloid-counterion complexes which are over- 
charged (charge-reversed) — the number of condensed coun- 
terions can actually be larger than is necessary to completely 
neutralize the colloidal charge H [H [H [H HI lH]. If 
this happens, the electrophoretic mobility of colloidal parti- 
cles will be reversed. While there are some theories which 
qualitatively account for this curious behavior, no fully pre- 
dictive approach is yet available |34]. In this paper we will 
use extensive Monte Carlo simulations to explore two aspects 
of this problem — the dependence of the minimum colloidal 
charge at which the reversal of the electrophoretic mobility 
takes place on: (one) the concentration of the multivalent z:\ 
electrolyte and (two) on the ionic size. 



II. THE MODEL AND SIMULATIONS 



The electrophoretic mobility, in general, is a complicated 
non-linear function of the electrokinetic C, potential ll36l[37ll . 
For small and large ionic strengths, however, the relation- 
ship between the two is linear and is given by the Smolu- 
chowski equation ll36i Issll . A change in the sign of the 
C, potential will, therefore, lead to the reversal of the elec- 
trophoretic mobility, which we will also associate with the 
overcharging (or the charge-reversal) of the colloidal particles. 
In principle, the overcharging (charge-reversal) and the rever- 
sal of the electrophoretic mobility are two distinct concepts - 
one is static and the other dynamical. In practice, however, 
the definition of charge reversal carries some ambiguity. The 
general trends, such as the behavior of the effective charge 
and of the electrophoretic mobility as a function of, say, the 
ionic size or electrolyte concentration are very similar in two 
cases ||39ll. We will, therefore, use the two concepts inter- 
changibly. For a fixed electrolyte concentration, the value of 
the colloidal charge at which the ^-potential vanishes will be 
designated as the critical colloidal charge. Our goal is to find 
the dependence of this charge on the ionic size and the elec- 
trolyte concentration. 

We consider a diluted aqueous mixture of colloidal parti- 
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cles inside a z:l electrolyte. The spherical colloidal particles 
have radius R surface charge —Zq, where q is the elemen- 
tary charge. For each colloidal particle there are Z monova- 
lent counterions. All the ions are modeled as hard spheres 
of radius a with the charge +zq (electrolyte counterion), —q 
(electrolyte coion) or +q (colloidal counterion) at their cen- 
ters. The solvent is treated as a continuum of dielectric con- 
stant e. The relative strength of the electrostatic interactions, 
as compared to the thermal energy, is measured by the ratio of 
the ionic radius to the Bjerrum length, = q^ /AiTekBT. 

As was argued in Ref. ll39ll . the value of the C poten- 
tial may be associated with the electrostatic potential at the 
effective shear plane removed from the colloidal surface by 
one ionic diameter |40]. Since the maximum of the electro- 
static potential also occurs at approximately the same posi- 
tion, the precise location of the shear plane does not influ- 
ence strongly the value of the (^-potential. In this respect, our 
approach is quite similar to the one adopted by Bjerrum for 
simple electrolytes W„ AZ] . The static potential at the ef- 
fective shear plane can then be calculated using the canoni- 
cal Monte Carlo (MC) simulations |39]. Recent simulations 
show the basic correctness of this picture for normal wetting 
surfaces iIstI f39ll . Working with ^-potential is also advanta- 
geous as compared to defining the effective charge in terms 
of condensed counterions located within a sheath surrounding 
the colloidal surface. Such definition carries a large degree 
of arbitrariness, since the condensed counterions will in turn 
drive a co-associations of coions. The effective charge will 
then be strongly sensitive to the precise value of the sheath 
width. This is not the case for (^-potential, which under the 
same conditions develops a maximum near the colloidal sur- 
face which diminishes its sensitivity to the precise location of 
the shear plane. Furthermore, since the (^-potential is calcu- 
lated by integrating the electric field over the whole space, it 
already takes into account the layering effect that hinders the 
geometrical definition of the effective charge. 

A colloidal particle is fixed at the center of a cubic simu- 
lation box of side length L and is surrounded by the coun- 
terions and coions, the number of which satisfies the over- 
all charge neutrality. We define C as the molar concentra- 
tion of the z-valent counterions derived from the dissociation 
of z: 1 (strong) electrolyte — assumed to be fully dissociated 
in an aqueous environment. The electrostatic interactions are 
computed using the Ewald summation method [43] with 518 
Fourier-space wave vectors and a real-space damping param- 
eter K ^ 5/L. 

Two types of MC moves were utilized — ion transfer to 
a completely new random position inside the simulation box, 
which is useful for low salt concentrations, and a small linear 
displacement for high salt concentration, in order to give the 
standard acceptance ratios for the Metropolis algorithm. The 
number of microions in each simulation was varied from ap- 
proximately 50 up to 3000 particles, depending on the molar 
salt concentration and the box length. Typical runs involved 
10^ Monte Carlo steps for equilibration and lO'* steps for pro- 
duction. After equilibration, the average number of counteri- 
ons and coions in concentric spherical shells of equal thick- 
ness around the colloid were accumulated in order to obtain 



the density profiles pi{r). The mean electrostatic potential at 
distance r from the colloidal particle is then calculated as 
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where E{r) is the electric field and P{r) is the integrated 
charge (in units of q) within a distance r from the center of 
the colloidal particle. 



P{r) = -Z 



4:TTr''^dr' , 



(2) 



where i refers to the type of the microion. Since the typical 
integrated charge rapidly decays to zero |39], the upper cutoff 
in Eq. ^ is taken to be L/2. Following the Ref. iH], the 
shear plane was located at one ionic diameter from the col- 
loidal surface, so that ( = where Rs — R + 2a. 



III. SCALING ANALYSIS 

The model presented in the previous section is quite com- 
plex, with a number of distinct length scales. To organize 
and interpret the data of the Monte Carlo simulations we 
shall, therefore, appeal to the dimensional and scaling anal- 
ysis. There are five basic length scales: R, L, a, Xb, and 
C^^/'^. Since we are interested in very dilute suspensions, 
L ^ oo. Although this limit can not be achieved in the simu- 
lations, our box size was always taken to be sufficiently large 
so that critical colloidal charge —Zcq did not have any explicit 
dependence on L. We are, therefore, left with four relevant 
length scale, so that for a fixed z : 1 electrolyte, Zc is a func- 
tion of only three dimensionless ratios 



(3) 



where /(x, y, z) is a scaling function. Furthermore, we note 
that when the Debye length Cr> = V \/^t^^b{z'^ + z)C, is 
sufficiently short, R/£,d ^1 — which is almost always the 
case near the isoelectric point — the curvature effects will 
be screened, and the critical colloidal charge must be propor- 
tional to the colloidal surface area. This means that 



Zr 



AttR^ 



Xb 



.XbC'/' 



(4) 



We conclude that for sufficiently large salt concentrations, 
the reversal of the electrophoretic mobility will take place 
when the modulus of the colloidal surface charge density, 
a = Zq/AirR^, is larger than the critical value ctc, which 
depends on the Bjerrum length, ionic radius, and the concen- 
tration of electrolyte through the scaling function g{x,y), 



(5) 



The similarity transformation, Eq. is particularly useful 
when one wants to obtain the critical surface charge den- 
sity for suspensions with large concentrations of electrolyte. 



In these cases, the direct MC simulations become extremely 
slow, due to large number of microions which must be used 
to simulate a dilute colloidal suspension in the L ^ oo limit. 
However, Eq. (|5]l tells us that this critical surface charge den- 
sity can also be obtained by simulating a much smaller system 
at a slightly lower temperature and with a somewhat larger 
microions. For example, suppose that we want to find the 
critical surface charge density of colloidal particles inside a 
dilute suspension at room temperature, = 7.2 A, con- 
taining 3:1 electrolyte at concentration C = IM, with ions 
of radius 2 A. Instead of doing the direct simulation of this 
system, we can simulate a "similar"system with say half the 
number of microions, C = 0.5 M, at a slightly lower temper- 
ature, \b = 7.2 X 2^/^ = 9.07 A and with ions of radius 
2 X 2^/^ = 2.5 A. From Eq. (|5]l, the critical charge of the 
original system (at concentration C — \ M) will be 2'^/'^ times 
the critical charge of the "similar" system. The latter simula- 
tions, however, are much easier to perform since the number 
of microions involved is much smaller 

In Refs. ll44ll and Isill it was argued that the critical surface 
charge density is determined by the work that must be per- 
formed to transfer a multivalent counterion from the bulk elec- 
trolyte to the colloidal surface. In particular, it was found that 
tTc ^ (A/j,)^, where A/i is the change in the ionic solvation 
free energy between the bulk and the colloidal surface jiill . 
In the limit of vanishing a, A/i diverges as 1/a. This di- 
vergence is a consequence of Bjerrurn pairing of oppositely 
charged ions in the bulk electrolyte IH |4ll 14511 . The criti- 
cal surface charge density, therefore, diverges as Cc ^ 
Thus, the scaling function with x = 0, g{{), y), should be con- 
stant for all values of y. For small, but finite values of x, 
we expect to see a deviation from this behavior in two limits: 
when the entropic effects begin to dominate the electrostatics 
C^/^Xb < 0.5, and when the hard-core repulsion begins to 
dominates everything, C^/"^a > 0.2. 



IV. RESULTS AND DISCUSSION 

In Fig. [U we show the (^-potential as a function of the col- 
loidal charge density for a suspension containing 3:1 elec- 
trolyte at a molar concentration of C — 0.1 M. The radius of 
the colloidal particle was fixed at i? = 30 A, while its charge 
Z was varied. Most of the simulations were performed with 
the boxes of side length L = 120 A, which is large enough 
to produce very small colloidal volume fractions, thus, min- 
imizing the influence of the periodicity on the ionic distribu- 
tion. However, for low salt concentrations we increased the 
box length to L = 150 A and 210 A, in order to increase the 
number of particles and to obtain a better statistics of ion dis- 
tribution around the coUoid. 

For weakly charged colloidal particles the increase (in 
modulus) of the surface charge density was accompanied 
by a uniform decline of the ( potential (C accompanied the 
colloidal charge and became more negative). However, when 
the colloidal charge became sufficiently large, counterion 
condensation became important and C increased as a function 
of the bare colloidal charge, becoming positive for sufficiently 
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FIG. 1: Zeta potential as a function of the colloidal surface charge 
density. The molar concentration of 3 : 1 electrolyte is C = . 1 M, the 
radius of the microions is a = 3 A, and the Bjerrum length is 7.2 A. 
The inset shows the region where the potential becomes reversed, 
( = 0, with a very good linear fit (solid line) to the simulation data 
(circles), from which the precise value of Uc is determined. 



strongly charged colloids, see Fig.[T] To accurately determine 
the critical colloidal charge density CTc at which C = 0, we 
used a linear interpolation of the simulation data, as shown in 
the inset of Fig. [T] 
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FIG. 2: Critical colloidal surface charge density as a function of the 
molar concentration of 3:1 electrolyte. Circles are for electrolyte 
with a = 2 A, As = 7.2 A; squares are for electrolyte with a = 3 A, 
A_B = 10.8 A. In both cases the ratio a/As ~ 0.278 is the same. 
Inset shows the data collapse when the concentrations and ac are 
properly scaled. 



We next studied the dependence of the critical surface 
charge density on the concentration of electrolyte. In Fig. 2, 
CTc is plotted as a function of C for two different electrolytes: 
a = 2A, As = 7.2 A; and a = SA, As = 10.8 A. Although 
clearly distinct, the two systems are "similar", since the ratio 
As/a is the same in both cases. Therefore, if OcO^jq is 
plotted as a function of XbC^^^ ( or as a function of CA^) 
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FIG. 3: Scaled critical colloidal surface charge density as a function 
of the scaled concentration of 3:1 electrolyte. The microions are of 
radius o = 3 A. 



the data for the two systems should collapse onto a single 
curve. This is precisely what is found, see the inset of Fig. 2. 

If the electrolyte systems are not connected by the 
similarity transformation, a data collapse is not expected. 
Nevertheless, as was discussed in the previous section, if 
a/As is small, g{x « 0,y) should be nearly constant. We 
expect, however, the deviation from this constancy to take 
place when the product xy becomes sufficiently large — the 
separation between the microions becomes compatible to the 
ionic size, or when the concentration becomes so small that 
the entropic effects dominate over the electrostatics. In Fig. 3 
we plot (TcO^/qasa function CAg, for various electrolytes of 
different Xb- 
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FIG. 4: Scaled critical coUoidal surface charge density as a function 
of the scaled concentration of 2: 1 electrolyte. The microions have 
radius o = 2 A. 



The figure clearly shows the inflection point of CTc as a 
function of the electrolyte concentration. Furthermore, in 
the limit of vanishing a/ As, the inflection point turns into 
plateau which extends up to the concentrations for which 
C^/^a « 0.2, consistent with the discussion presented earher. 
In the limit of vanishing ionic size, the plateau region extends 
indefinitely. 



V. CONCLUSIONS 

Using extensive Monte Carlo simulations we have studied 
the dependence of the minimal colloidal charge at which the 
reversal of the electrophoretic mobility first takes place, on 
the concentration of 3:1 electrolyte and on the ionic size. The 
critical surface charge density CTc was found to exhibit an in- 
flection point as a function of the electrolyte concentration. In 
the limit of small a/Xs, the inflection point becomes a flat 
plateau, extending from the lower concentration Ci to the up- 
per concentration C„. The value of the lower bound is de- 
limited by the distances at which the entropic effects begin 

1/3 

to dominate over the electrostatics, C; Xb ~ 0.5, while the 
value of the upper bound is determined by the distances at 
which the hard-core repulsion begins to dominate over ev- 

1 /3 

erything C« a « 0.2. In the interval [Ci,Cu\ the surface 
charge density is found to be CTc ~ 0.013g/a^, independent 
of the electrolyte concentration. The plateau disappears, turns 
into a simple inflection point with vanishing second deriva- 
tive, when Xb/o, < 2.5. The same behavior was found to 
occur for electrolytes of other valences. For example, in Fig. 
4 we plot the critical surface charge density as a function the 
electrolyte concentration for 2: 1 electrolyte. Once again the 
inflection point and the formation of the plateau are evident. 
However, in this case, to compensate for the weaker electro- 
static interactions between the coions and the counterions, the 
value of Ab must be significantly lower for the plateau to ap- 
pear clearly. What is surprising, however, is that the scaled 
surface charge density (JcO^/q, appears not to depend on the 
valence of the electrolyte — or depend only very weakly — 
in the plateau region. In the case of 2: 1 electrolyte we find 
fliat acO^/q « 0.015 as compared to the 0.013 for the 3:1 
electrolyte. It will be interesting to see if this curious be- 
havior persists for other values of z. At the moment, there 
is no theory which can quantitatively account for these curi- 
ous findings. We hope that the present study will provide a 
simulational benchmark against which the future theoretical 
predictions can be tested. 
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